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ABSTRACT 

The  Schwinger  variational  method  is  formulated  for 
three-body  collisions.  When  exchange  effects  can  be  excluded, 
the  generalization  from  two  to  three  body  problems  is  straight- 
forward. Exchange  effects,  however,  introduce  difficulties  in 
the  latter  problem  lAiich  do  not  occur  in  the  former.  Despite 
this,  the  method  can  be  extended  to  give  a  stationary  expres- 
sion for  the  exchange  scattered  amplitude.  As  a  result  of 
this  formulation,  one  gains  some  insight  into  a  superior 
method  for  treating  three  body  problems  with  exchange  in 
the  Born  approximation. 
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I.   Introduction 

A  controlled  eaq)erimental  investigation  of  the  physical  processes 
in  the  ionosphere  presents  enormous  difficulties  inasmuch  as  the  conditions 
which  exist  there  cannot  he  reproduced  terrestrially.  The  hest  one  could 
hope  for  would  be  to  isolate  some  physical  processes  to  he  studied  on  earth. 
In  order  to  apply  the  results  of  such  experiments  to  ionospheric  data,  it  is 
important  to  know  on  which  physical  parameters  the  experimentally  measTired 
quantities  depend,  and  what  this  dependence  is.  A  considerable  amount  of 
information  could  he  obtained  from  a  theoretical  investigation  of  these 
fundamental  processes.   If  one  finds  that  the  results  of  terrestrial  ex- 
periments can  be  predicted  theoretically,  one  can  with  some  assurance  extra- 
polate to  the  conditions  existing  in  the  ionosphere.  Thus  theoretical  in- 
vestigations are  of  some  importance  both  as  a  guide  to  ea5)erimental  work  and 
for  the  prediction  of  ionospheric  effects. 

The  difficulty  in  the  theoretical  approach,  however,  lies  in  the 
fact  that  even  the  simplest  problems  involve  many  particles.  Just  as  in 
classical  mechanics,  quantiun  theory  does  not.  In  general,  permit  any  rigor- 
ous solution  for  this  type  of  problem.  A  program  for  the  theoretical  in- 
vestigation of  the  physical  processes  in  the  ionosphere  must  necessarily 
begin  with  a  mathematical  investigation  of  suitable  approximation  methods. 

In  this  report  such  a  program  is  begun  in  connection  with  scatter- 
ing problems.   This  does  not  constitute  by  any  means  the  first  attempt  for 
the  solution  of  problems  of  this  type.  Mott  and  Massey  and  their  collabora- 
tors have  solved  many-particle  collision  problems  of  all  sorts  using  princi- 
pally the  Bom  approximation.  This  method  assumes  that  the  incoming  wave  is 
distorted  very  little  by  the  scattering  center. 

Because  of  the  low  temperature  of  the  ionosphere  the  particles 
there  possess  energies  of  only  a  fraction  of  an  electron-volt.   Under  these 
circTJmstancee,  the  scattered  particle  spends  so  much  time  in  the  neighbor- 
hood of  the  scatterer  that  the  approximation  of  small  distortion  of  the  In- 
cident particle  is  no  longer  valid. 


1.  N.  ?.  Mott  and  H.S.W.  Massey,  Theory  of  Atomic  Collisions,  Ch.  8,  Oxford 
University  Press. 
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An  additional  difficulty  arising  in  three  "body  collisions  is  in  the 
possihility  that  if  the  incident  particle  is  identical  with  one  of  the  particles 
in  the  scatterer,  exchange  collisions  may  take  place;   i.e.,  particle  1  may 
come  in  and  he  captured  and  particle  2  emerge.  According  to  the  Pauli  prin- 
ciple, this  type  of  scattering  cannot  he  distinguished  from  direct  scattering 
in  which  the  incident  particle  1  is  itself  scattered.  Exchange  collisions 
may  "be  of  considerahle  importance  for  low  energy  encounters;  and  yet  the 
theory  of  exchange  collisions  in  any  known  approximation  is  not  "based  on 

sound  theoretical  foundation.  Calculations  made  according  to  existing  pro- 

2 
cedares  seem  to  he  in  serious  disagreement  with  experiment  . 

Only  recently  have  attempts  "been  made  to  find  other  more  suitable 
^proximations.  A  variational  method  of  Hulthen  has  "been  generalized  to  in- 
clude three-hody  collisions  .  Elastic  and  inelastic  processes  have  also  "been 

H  5  6 
treated  extensively.'  • 

The  success  of  the  application  of  the  Schwinger  variational  proced- 

7  8 
ure  to  two-body  collisions  *  has  made  its  generalization  to  three-body  problems 

seem  worthwhile.  Another  variational  principle  formulated  in  terms  of  G-reen's 
functions,  while  more  difficult  to  apply,  gives  more  accurate  results.  This 
report  attempts  to  formalate  such  a  variations^  principle  which  would  be  use- 
ful in  discussing  the  scattering  of  electrons  by  hydrogen  atoms  both  when  ex- 
change scattering  is  included  and  when  it  is  not  included.   Some  remarks  will 
be  made  about  the  difficulties  encountered  when  identical  particles  are  in- 
volved and  about  the  general  usefulness  of  the  Schwinger  variational  method 
in  problems  of  this  sort. 


2.  Bates,  Hassey,  Fundaminsky  and  Leech,  Phil.  Trans,  of  Roy.  Soc.  A, 2^3,93  (195O) 

3.  ',/.  Kohn,  Phys.  Rev.  JU,  I793  (191^8). 

U.  Su  Shu  Huang  -  Phys.  Rev.  76,  ^77  (19I+9) 

5.  H.  S.  ¥.  Massey  and  B.  L.  Moiseiwitsch,  Proc.  Roy.  Soc.  A205,  U83  (I95I). 

6.  M.  Verde,  Helv.  Phys.  Acta  XXII,  339  (19^9);  A  Troesch  and  M.  Verde,  Helv. 
Phys.  Acta  XXIV  39  (I95I). 

7.  J.  Schwinger,  Lecture  Notes  on  Nuclear  Physics,  Harvard  University. 
J.  Schwinger,  B.  A.  Lippman,  Phys.  Rev.  79,  kSl   (I95O). 

g.  Blatt  and  Jackson,  Rev.  of  Mod.  Phys.  22,  77  (I95O). 
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II.  Exchange  Excluded 


\'!e   consider  in  this  section  a  collision  of  particle  1  incident  on  a 
hound  system  of  particles  2  and  3.  All  three  particles  are  different,  so  the 
complications  arising  from  the  Pauli  principle  can  he  neglected.   The  mass  of 
particle  3  is  taken  to  he  infinite.  This  assumption  is  suitahle  for  a  dis- 
cussion of  the  scattering  of  electrons  hy  hydrogen  atoms  with  exchange  effects 
neglected.   It  would  not  he  suitahle  for  the  description  of  the  analogous  nuclear 
prohlem,  i.e.,  the  scattering  of  neutrons  hy  deuterons.  For  this  latter  prohlem 
the  generalization  of  the  work  in  Section  H  is  straightforward,  hut  the  material 
in  Section  m could  not  he  generalized  so  easily. 

In  all  scattering  prohlems  the  calculated  quantity  which  is  confronted 
vfith  experiment  is  the  differential  scattering  cross-section.  This  is  defined  as 
the  current  of  scattered  particles  crossing  unit  area  at  infinity  at  some  par- 
ticular angle,  divided  hy  the  current  of  particles  per  unit  area  in  the  incident 
heam.   This  miraher  can  he  found  hy  considering  the  solution  of  the  Schrodinger 
equation,  which  is  the  following  for  this  prohlem: 

In  this  equation  r.  and  r  are  the  position  vectors  of  particles  1  aT>H  2  respect- 
ively. The  position  of  particle  3  is  taken  as  the  origin  of  coordinates.  The  Vs 
are  the  potentials  hetween  the  ^articles,  and  we  shall  assiime  that  the  forces  axe 
centred;  E  is  the  total  energy  of  the  system,  i.e.,  the  sum  of  the  kinetic  energy 
of  the  incident  particle  and  the  energy  of  the  particle  in  its  initial  state  (which 
is  negative). 

The  asymptotic  form  of  the  solution  of  (2.1)  should  represent  an  in- 
coming wave  plus  outgoing  spherical  waves.  If  we  assume  that  particle  1  is  incom- 
ing  along  the  negative  z  eixis,  the  asymptotic  form  of  the  solution  is 

ik(a.r,  )  ^K\^^\ 

L    ^(r^,r2)  =  e  ""  ^"^     ^  0^(r2)+ H  V^'S^^^^  TTT^  K^'^P^ '  ^^.2) 

r.-»oo  n  I  l| 


♦In  what  follows  all  summations  include  integrations  over  continuum  states. 


.  l^- 


Here  0  are  the  wave  functiona  of  the  stationary  state  of  the  potential  Vg-rCr-), 


2_  2^ 


j^=  _  (E-E^).  and 


/^s 


Is  a  unit  vector  along  the  z  axis,  the  direction  of  the  Incoming  electron. 


£  is  the  energy  of  the  n   state  of  system  2  and  3. 

If  the  asymptotic  form  of  the  solution  is  given  ■by(2.l),  the  scattering 
cross-section  for  excitation  of  the  hound  system  to  the  n   state  can  he  shown 
to  he  k     2 

<r„=  E^ig  .  (2.3) 

0 

In  order  to  derive  a  Schwinger  variational  method,  it  is  necessary  to 
replace  equation  (2.1)  hy  an  Integral  equation  which  satisfies  the  houndary  con- 
dition (2.2).  To  achieve  this  result  we  transpose  "^3^3(1"^)  ^^'^   ^12^'^1»'2^  in (2.1) 
to  get 

r[^V=  -^V^].CE-7„)W=(T„.7i,)V=TV/.        C2.U) 

where  we  have  written  V  for  (V^,  +  '^Ig^* 

The  integral  equation  replacing  (2.1)  and  satisfying  (2.2)  is  given  "by 

(2.5) 
¥°   ^^V^2^   =  «^  ^K   ^/'z-'l^  ^o^'2^  +/dr|,dr.  ^(r^.r^.  r^.^^^^il^^K^^^''^^' 

In  equation  (2.5)  o„ 

I   I    2m^   ,  I   I, 
v{t^,   Tg)  =  — ^  V(rj^,  T^   . 
a. 

y/^     is  a  solution  of  (2.1)  which  asymptotically  represents  an  incoming  plane 

p 

wave  of  particle  1  with  energy  k  and  a  hound  state  of  particle  2  with  energy  E^. 


•   I. 


G(r^,r  ;  r   r.)  Is  the  Green's  function  of  the  operator 

and  satisfies  the  differential  equation 


{\[k^l*h^^ly-^^^} 
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N 


*  IC  -  V23 


.  a{T^,T^;r^,T^)  =  +8  ir^^-r^^)   6(r2-r2')  .  (2.6) 


The  (Jreen's  function  can  "be  fo-und  "by  the  method  of  separation  of  vari- 
ahles  to  "be  explicitly 

G  =  -  T"    0     (r    )0(r)         ■      ,    ^^      ,^ 

*  exp  ik   Itt  -r.    |  ^2.?) 


This  result  follows  from  the  fact  that 

H^n  ^'2^  0^(r2)  =  8(r2-r2  ),  (2.g) 

exp  ikjr^-r^^l 
and  that     Im  I  r   -  r  M        ^^  *^®  Green's  function  for  the  operator  (Vf^  +  k    ). 

In  the  limit  r^**  00   ,  (2. 7)  "becomes 


r- ->oo 


exp  ik   I  r.  I  1*1 

^-'^K^-2^     pTp—     -^^-P-  V/"-    -l)5*n('2^-      (2.9) 

Inserting  (2.9)  into  (2. 5^and  comparing  with(2.2)we  have 

^"^  ^n  =  -y*exp-ik^  (^-  r^^")  0^'  (rj)  v(rj_'  r^)^^^  (r^',  r^)  dr^^  drj  .  (2.10) 

for  the  scattered  amplitude.  The  cross  sections  can  he  determined  from  (2.10)  accord- 
ing to  equation  (2.3). 

In  order  to  formulate  a  variational  principle  for  f''  ,  it  is  necessary  to 
write  down  an  integral  ecjuation  adjoint  to  (2.5)  which  satisfies  (2.U).  This  one 
is 

^^(ri.r^)  =  exp  ik^(^-  r^)  ^^(rj)*^ G^r^.r^)  (r^',r2)v(r^'.r2)^^r^',r2)  dr^^"  dr^'.  iZM) 
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The  asymptotic  form  of  this  solution  consists  of  concentrating  waves.     By  multiply- 
ing the  complex  conjugate  of   (2.1l)  ty  vj^     and  integrating,    one  sees  that 

Utt  f°  -fyp-l  r^)  v(r;  rj)^  (r^.r^)  dr^  dr^  (2.12) 

7/ f-  ^"l'"2^^("l'2^  '^^^1''2   '    ^1  '2^  -^'l''^2%^^'^l''^2)<^1^2^1^2    ' 
By  treating  Equation  (2.5)   in  an  analogous  way  we  have 

im  f^  ^y^exp  ik^C/^-r,)  ^.CrJ)  v(r;.r2')  jj^^  (r^'^r^')  dr,'  dxj  =  hv  f°.        (2.13) 

q 
Vith  the  reciprocity  relation  (2,13),    one  can  form  the  expression-^ 


imf' 

n 


(2.1U) 


-Jjf>l''2^  v(r^.r2)G(r,.   r^^T^,r^>r{T^,T^):^^{r^,v^)^^^^  6x^  dr^     . 

~^^~~~~  *n  /"  0  ' 

yexp  iko(/^2.  Tj^')  ^^(rj)  v(r^',r2)  ^^  (r^^^.r^')  dr^^  dr2yexp-ik^yx.r^')0*(r2)j^^(^'i^  )d^^' 

which  is  stationary  with  respect  to  arbitrary  variations  oty     and  y/     ,  Thus  if 

*n  /%  f 

trial  functions  were  substituted  for  W     and  ^  on  the  right  hand  side  of  (2,  lU) 

which  differed  from  the  true  solution  by  ^  ,  the  scattered  amplitude  f  ,  would  be 
correct  to  order  £    ,  obtained  from  (2.lU). 


9.  N.  Marcuvitz  -  Sec.  III.  D.  "Recent  Developments  in  the  Theory  of  Wave  Propaga- 
tion", New  York  University  Institute  for  Mathematics  and  Mechanics,  19^9. 
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III.  Exchange  Included 

If  particles  1  and  2  are  Identical  the  results  of  Section  II  are 
not  complete.  Under  these  circumstances  we  have,  as  indicated  in  the  Intro- 
duction, the  possihility  of  exchange  scattering,  namely  the  possihility  that 
particle  2  and  1  will  change  roles  and  the  former  will  emerge  when  the  latter 
is  incident. 

This  exchange  scattering  will  now  affect  the  value  for  the  total 
cross-section  for  the  following  reason.  The  Paoli  principle  imposes  certain 
symmetry  conditions  on  the  wave  functions.  The  complete  solution  (space  and 
spin  parts)  must  he  ant i symmetrical  for  Jermi  particles  and  symmetrical  for 
Bose  particles.  Since  the  applications  we  envisage  are  concerned  with  electron 
scattering,  we  shall  restrict  ourselves  to  the  former  case,  in  which  the  space 
part  must  he  anti-symmetriceil  if  the  spins  are  parallel,  and  the  space  part 
mast  he  symmetrical  if  the  spins  are  anti-parallel.  It  can  he  shown   that 
the  quantities  which  determine  the  total  differential  scattering  cross-section 
are 

and  (3.1) 

-  1^      |2 
n   '  n   Ti' 

f  is  the  ordinary  scattering  cross-section  as  determined  in  Section  II  ,  and 

g^  is  the  exchsinge  scattering  cross-section  \^ich  we  have  discussed  ahove. 

The  total  scattering  cross-section  is  then  given  hy  a  suitahle  weighted  mean 
of  the  quantities  in  (3.1),  namely 

k 


^n=  iT^^  ^*  ^  -n^-  (3.2) 

0 


It  is  possible,  therefore,  to  satisfy  the  Pauli  exclusion  prin- 
ciple without  explicitly  symmetrizing  the  result.  This  requires,  however, 
that  we  know  of  the  exchange  scattering  amplitude,  g^  . 

10.  Mott  and  Massey,  loc.  clt.  Chap.  8. 
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The  forcpalation  of  the  protlem  given  in  Section  II  is  not  con- 
venient  for  deecrilDing  the  possibility  that  particles  1  and  2  exchsmge  their 
roles.      In  addition,   Kohn       has  pointed  out   that   it  is  not  certain  whether 
any  iterative  procedure  Tsased  on  equation  (2.5)  could  descri"be  a  "bound  state 
in  particle  1  and  a  continuum  state  in  particle  2.     Both  difficulties  can  "be 
ciroiravented  fey  treating  "both  particles  in  a  symmetric  way.      We  rewrite  equa- 
tion (2.1) 


v;here 


,2  .„2 


V^  +y^+  e  -  Vi,(ri)  - '^^-h^^'z^ 


2m 
2m 


'23' 


f=  V  (r^.r^)]^ 


(3.3) 


r(r     r   )     =     —    V 
■"■  ^  ^  ^2       12 


'13 


2mT 
i.  V 


2m 
23  ^2  '23 

The  integral  equation  corresponding  to    (5.3)  is 


)^=  0^(r^)0^(r2)  ^J(i(v^,T^^vlv'^h^^{Tl,Tl)l|Ji.Tly^)  dr^'  drj  (3.U) 


I      I 


where  G(r- ,rp;r- ,r^)   is  the  Green's  function  of  the   operator 

L  = V^  +V^+e-v       -V 
^         1       ^2  13         23 

and  0.  (r  ),  ^  (r)  are  solutions  of  the  equations 

^'7'*«o--23^95,  =  0 

(v2+(e-e^)-  v^^,)  0^  =  0; 


(3.5) 


(3.6) 


e  =  ground  state  energy  of  the  system  consisting  of  particles  2  and  3. 


11.  W.  Kohn  -  private  communication. 
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In  order  that  the  second  part  of  (3.U)  give  the  scattering  cross- 
section  directly,  ^  at  infinity  should  represent  a  wave  coming  in  along  the 

p 
direction  of  particle  2,  with  energy  k  =  e-e  . 

0       0 

The  Green's   function  has  two  representations  which  again  can  easily 
"be  found  "by  separation  of  variables  using  the  relation  (2.8): 

GCr^.r^;   rj^^,r2)  =  -i:    0^(r^)0;(r^')     g^^       (r^.r^)  (3.7a) 


=  -  H    {^n^^2^  K  ^^^2^  e^     (r^.r^)  .  (3.7I3) 

n 

If  v,^  falls  off  faster  than     —  (or  — -  )    ,    the  asymptotic  form  of  (3.7a) 

12  ^      2 

with  r^  ->oo  is 

and  the  limit  as  r^->oo  can  he  fotind  from  (3.7"b): 

^^nl^'l' 

In  (3.8a)  and  (3.-ffb)  F^^  represents  the  solution  of  the  homogeneous  equation 

n 

(3.6)  which  is  an  incoming  wave  at  infinity  in  the  direction  TT  -0or  (77  -  i>) 

with  energy  e-e^.  0  and  ^  are  the  angles  between  r^   and  t       and  r.  and  r. 

Using  equations  (3.8a)  and  (3.ffb),  we  can  now  find  an  expression 
for  f°  and  g^  from  (3.1^){ 


/7 


,       ....  (3.9) 


♦The  Coulomb  case  will  be  discussed  in  Sec.    IV. 
12.  Mott  and  Massey,   loc.cit.   Chap.^. 
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and  ^^n'^1 


) 


^2^^  ^  1^1  ^  (3.10) 

W//      »        i>  I  I 


thus 


and 

4  =  -  ^fl>\4K  ^^1  ^-^^  ^12^^1''2V^V^2)   -^l   ^2   •  ^5.12) 

J  n 

It   should  "be  noted  here  that  the  elastic  scattered  amplitude  is  not  given  com- 
pletely "by  f°   calculated    according  to    (3.11 ).     The  outgoing  spherical  part  of 

jj.  (r,)  0   (Tp)  should  "be  added  to  this  value. 

As  in  Section- II ,  we  can  now  write  down  stationary  expressions  for 
l/kT\f°  and  l/iffTg°  .      We  note  that   the   equation  adjoint  to    (3.^)  for  the  calcula- 


tion of  f°  is 
n 


^^(r^.r^)  =  \  (r2,TT-0)jS^(r^)  +  /o'Cr^.r^^rj^.r^)  Y^^{T^,T^)')({T^,r^)  dr^  drj   .(3.13) 
n  ^ 


The  a^djoint  equation  for  the  calculation  of  g     is 


^(r^.r^)  =  F*  iT^,v-^)p^{T^)  *fG\T^,T^iT^,r^)r^^iT^,T^)  Mr^,T^)  dx^  drj  ,     (3.11+) 

n  *' 

The  reciprocity  relations  corresponding  to  (2.13)   can  again  "be  proved  in  the   same 
way,  and  we  have  for  f     and  g^  following  expressions  which  are  stationary  with 
respect   to  arhitrary  variations  of  ^,  X,  and  ](.  i 

It      II  n      It  I  I 

1  "  v^2^'^l''^2^^^'^l'^2^        '    "^1     ^2       <^i.^2 


l+nf 
n 


a      /s*i(^l)i^o^'"2^''l2^'"l'''2^'^*^'^l«''2W  ^2F]L^1'^-^^'^*^''2\2^V2^1^^1'''2' 


I        I 

dXj^  dTg 
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and  /'    *,     '      '\         /'      '\  //'      '\       /'/'*/'      'n        /     '      'n    /     '      '      "      " 


(3.16) 


t       I 


^(r^,r2)  dr^  dr^ 

Having  evaluated  f     and  s     we  can  find  the  cross-section,   including 
exchange,    for  spin  ^particles  "by  reference  to   (3.1)  and  (3.2). 


IV.   The  Coulomti  Case. 

If  the  potentials  v,,  and  Vp-  go  to   zero  as  and  respectively, 

as  would  "be  the  case  if  we  were  discussing  the  scattering  of  electrons  "by  hydrogen, 

then  some  of  the  statements  made  in  the  preceding  section  are  not  precisely  correct. 

iklr^l  ik|rJ-i]og  k|r3_| 


6  6 

We  have  the  usual  difficulty  that  the  outgoing  is  not  hut  — 


and  there  are  similar  difficulties  with  the  incoming  wave.''  We  may  carry  through  ex- 
actly the  same  program  as  we  did  in  Section  III  after  investigating  the  asymptotic 
form  of  the  hydrogenic  Green's  function  and  after  the  correct  interpretation  of  the 
scattering  cross-section  has  "been  made  for  this  case.  The  formulae  derived  are  the 
same,  provided  we  choose  the  0,  and  F  functions  to  represent  the  incoming  waves 
corresponding  to  the  Coulomh  case.  While  no  formal  difficulties  exist,  the  ques- 
tion of  the  convgrgence  of  an  iterative  procedure  for  the  Coulomh  field  is  still 
unsettled.   This  matter  is  rather  conplicated  end  v;ill  not  he  discussed  in  this 
report. 
V.  Evaluation  of  the  Results 

Having  carried  through  the  formulation  of  the  Schwinger  variational  method 
for  three  "body  collisions,  it  is  fair  to  ask  v;hat  are  the  prospects  for  utilizing  it 
in  the  coiapiitation  of  cross-sections  which  are  of  interest. 

13.  Mott  and  Massey  loc.cit,  p.  HS. 
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At  present  it  seems  that  it  is  necessaTy  to  treat  the  two  par- 
ticles in  a  symmetrical  way;     however,   this  requires  the  use  of  a  Green's  func- 
tion so  complicated  that  the  difficjlties  involved  in  its  use  would  he  almost 


prohihitive.     There  are,  however,   some  encouraging  aspects.     If  the  conjecture 

3ni 
2 


of  Kohn   is  correct,  it  may  "be  v/orthwhile  to  examine  the  symmetrical  treatment 


in  the  Born  approximation.  This  may  remove  some  of  the  previous  difficulties 
encountered  in  the  CEuLculatioas  including  exchange.   As  a  final  heartening 

note  if  one  excludes  exchange  effects  and  follows  the  formulation  of  Section 

lli 
II,  the  integrals  involved  are  tractable  '.  In  fact,  using  a  Born  trial  field, 

the  integrals  can  he  explicitly  evaluated  except  for  some  minor  approximations, 
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lU.  An  application  of  this  sort  vill  "be  the  suhject  of  a  future  report  hy 
H.  Boyet  and  S.  Borowits. 
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